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Abstract 

Supercritical  conditions  designate  a  situation  where  the  working  fluid  pressure 
is  above  the  critical  point.  Among  these  conditions,  it  is  interesting  to  identify  a 
transcritical  range  which  corresponds  to  cases  where  the  pressure  is  above  the  crit¬ 
ical  point,  but  the  injection  temperature  is  below  the  critical  value.  This  situation 
is  of  special  interest  because  it  raises  fundamental  issues  which  have  technological 
relevance  in  the  analysis  of  flows  in  liquid  rocket  engines.  This  situation  is  here  en¬ 
visaged  by  analyzing  the  behavior  of  a  nitrogen  shear  coaxial  jet  comprising  an  inner 
stream  injected  at  temperatures  close  to  the  critical  temperature  and  a  coaxial  flow 
at  a  higher  temperature.  Experiments  are  carried  out  both  in  the  absence  of  exter¬ 
nal  modulation  and  by  imposing  a  large  amplitude  transverse  acoustic  field.  Real 
gas  Large  Eddy  Simulations  are  performed  for  selected  experiments.  The  combina¬ 
tion  of  experiments  and  calculations  is  used  to  evaluate  effects  of  injector  geometry 
and  operating  parameters.  Calculations  retrieve  what  is  observed  experimentally 
when  the  momentum  flux  ratio  of  the  outer  to  the  inner  stream  J  =  (pew^)  /  (piwf) 
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is  varied.  Results  exhibit  the  change  in  flow  structure  and  the  development  of  a  re¬ 
circulation  region  when  this  parameter  exceeds  a  critical  value.  The  instantaneous 
flow  patterns  for  different  momentum  flux  ratios  are  used  in  a  second  stage  to  char¬ 
acterize  the  dynamical  behavior  of  the  flow  in  terms  of  power  spectral  density  of 
velocity  and  density  fluctuations.  Results  obtained  under  acoustic  modulation  pro¬ 
vide  insight  on  mixing  enhancement  of  coaxial  streams  with  a  view  of  its  possible 
consequences  in  high  frequency  combustion  instabilities.  It  is  shown  in  particular 
that  the  presence  of  strong  acoustic  modulations  notably  reduces  the  high  density 
jet  core  length,  indicating  an  increased  mixing  efficiency.  This  behavior  is  more 
pronounced  when  the  jet  is  placed  at  the  location  of  maximum  transverse  velocity 
fluctuations. 


Nomenclature 

a  Thermal  expansion  coefficient 

f3  Isothermal  compressibility  coefficient 
A  Thermal  conductivity 

H  Dynamic  viscosity 

fj1  Turbulent  dynamic  viscosity 

z/  Turbulent  kinematic  viscosity 

uj  Acentric  factor 

p  Density 

p*  Normalized  density  (p  -  Poo)/(pc  -  Poo) 
Poo  Far  field  density 
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Poo  Reservoir  initial  density 

pc  Time  averaged  centerline  density 
pe  Outer  jet  injection  density 

Pi  Inner  jet  injection  density 

r  Viscous-stress  tensor 

r*  SGS  stress  tensor 

a,  b ,  c  Coefficients  for  the  Peng-Robinson  equation  of  state 

cp  Constant  pressure  specific  heat 

cv  Constant  volume  specific  heat 

Cw  WALE  model  constant 

di  Inner  jet  diameter 

E  Total  energy 

es  Internal  energy 

Ep  Spreading  rate  of  the  density  profile 

he  annular  duct  thickness 

hs  Enthalpy 

J  Momentum  flux  ratio 

li  Inner  lip  thickness 

Lp  ’Half  Width  at  Half  Maximum’  of  the  density  profile 

/ o.i  t 

Lp  Xp  XQ 
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/0.5  ™0.5  ™ 

V  Xp  ~  X0 

/0.99  ™0.99  ™ 

Lp  Xp  ~  X0 

0  %u  l u 

p  Pressure 

Pc  Critical  pressure 

R  Ideal  gas  constant 

r  Specific  gas  constant 

S  Density  ratio 

S  Symmetric  part  of  the  velocity-gradient  tensor 

s  Symmetric  traceless  part  of  the  square  of  the  velocity-gradient  tensor 
T  Temperature 

Tc  Critical  temperature 

Tr  Reduced  temperature 

Reservoir  initial  temperature 
Taver  Duration  of  the  averaging  procedure 

Te  Outer  jet  injection  temperature 

Tj  Inner  jet  injection  temperature 

U  Velocity  ratio 

u  Axial  velocity 

u*  =  u/ue  Normalized  velocity 
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Uoo  Reservoir  initial  velocity 

ue  Outer  jet  injection  velocity 

Ui  Inner  jet  injection  velocity 

x^1  Axial  centerline  position  where  p  =  0.1(p;  —  pe)  +  pe 

Xp'5  Axial  centerline  position  where  p  =  0.5 (pi  —  pe)  +  pe 

x°p  99  Axial  centerline  position  where  p/ pi  =  0.99 

xu  Axial  centerline  position  where  axial  velocity  is  minimum 

qf  SGS  heat  flux  vector 

q  Heat  flux 

v  Velocity  vector 

w  Vector  of  conservative  variables 

PH  Turbulent  Prandtl  number 

Re  Reynolds  number 

SO  (fue)/(dt  +  lt) 

Ste  ( fue)/he 

St*  ( fui)/di 

Stz  (fue)/h 
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1  Introduction 


Flows  formed  by  coaxial  injectors  raise  many  fundamental  issues.  Such  flows  are  also 
of  practical  relevance  for  high  performance  cryogenic  liquid  rocket  engines  (LREs).  It 
is  important  to  understand  the  basic  processes  taking  place  in  these  flows  to  improve 
the  design  methodology  for  future  cryogenic  LREs.  Another  potential  application  of 
research  reported  in  this  article  lies  in  the  area  of  combustion  instability.  It  is  known 
that  mechanisms  driving  combustion  oscillations  result  from  the  coupling  between  the 
injector  flow,  the  combustion  process  and  the  chamber  acoustic  modes.  The  manner  in 
which  the  inner  and  outer  streams  originating  from  a  coaxial  element  interact  with  each 
other  and  with  their  surroundings  may  affect  the  temporal  and  spatial  release  of  chemical 
energy  within  the  chamber  leading  to  undesired  oscillations  associated  with  a  resonant 
acoustic  motion.  In  most  LREs,  performance  has  been  enhanced  by  augmenting  the 
chamber  pressure.  Current  LOx/H2  engines  like  the  SSME,  Vulcain2  and  RS-68  operate 
at  values  exceeding  the  critical  pressure  of  their  individual  propellants.  The  inner  jet, 
oxygen,  is  generally  injected  in  a  transcritical  or  liquid-like  state  (be.  at  a  pressure 
higher  than  the  critical  value  but  at  a  temperature  lower  than  the  critical  temperature) 
and  it  is  surrounded  by  a  higher-velocity  hydrogen  stream  at  supercritical  temperature 
[72],  The  inner  fluid  stream  injected  in  a  transcritical  form  evolves  to  a  supercritical 
state  as  its  temperature  rises  because  of  mixing  and  possibly  combustion.  Under  such 
thermodynamic  conditions,  the  system  exhibits  particular  features  that  differ  from  those 
of  a  two-phase  flow  or  of  a  pure  gaseous  injection.  The  objective  of  the  present  work 
is  to  compare  experiments  and  simulations  in  which  shear  coaxial  jet  configurations  are 
established  at  pressures  above  the  critical  point  of  the  injected  fluid,  but  at  transcritical 
temperatures  without  and  with  acoustic  excitation.  For  this  study,  nitrogen  (N2)  is  the 
only  fluid  used  in  the  inner,  outer  and  chamber  flows. 

At  this  point  it  is  worth  reviewing  previous  work  dealing  with  supercritical  injection. 
Starting  with  single  round-jets,  experiments  have  been  carried  out  at  AFRL  (Air  Force 
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Research  Laboratory)  [9]  and  DLR  (Deutsches  Zentrum  fur  Luft-  und  Raumfahrt)  [29,  41] 
(see  also  Oschwald  et  al  [43]  for  a  review  up  to  2006).  Some  recent  investigations  are 
reported  by  Segal  and  Polikov  [61]  and  Roy  and  Segal  [56].  These  studies  indicate  that 
the  liquid  break-up  and  atomization  mechanism  which  prevail  at  subcritical  pressure  are 
no  longer  observable  because  surface  tension  and  latent  heat  of  vaporization  vanish  at 
supercritical  pressures  [47].  The  jet  of  fluid  then  dissolves  in  the  ambient  gas,  with  no 
evidence  of  droplets.  The  flow  features  “comb-like”  structures  at  the  edge  of  the  dense 
stream,  a  type  of  pattern  which  is  not  observed  at  subcritical  pressures  [9] .  The  geometry 
of  the  flow  is  reminiscent  of  that  of  a  variable  density  turbulent  gas  stream  [9].  This 
behavior  has  been  confirmed  by  experimental  measurements  of  visual  spreading  rates  in 
the  initial  region  which  are  found  to  be  consistent  with  theoretical  growth  rate  expressions 
for  incompressible  variable-density  turbulent  mixing  layers.  Quantitative  measurements 
of  density  [4,  41,  29],  spreading  rate  and  axial  decay  coefficients  obtained  by  DLR  [43,  40] 
using  spontaneous  Raman  scattering  also  convey  the  same  type  of  information. 

Coaxial  injection  is  investigated  by  several  researchers  at  AFR.L  [14,  13,  25,  24,  55,  54] 
who  examine  the  single-species  case,  and  by  Mayer  and  Smith  [30]  who  examine  a  stream 
of  nitrogen  injected  at  low  temperature  (high  density)  surrounded  by  a  high  velocity 
stream  of  moderate  temperature  helium.  Quantitative  measurement  of  species  density 
are  reported  by  Oschwald  et  al  [42]  for  a  coaxial  nitrogen  /  hydrogen  injection,  using 
spontaneous  Raman  Scattering.  When  compared  with  the  round  jet  configuration,  a 
strong  reduction  of  the  jet  density  potential  core  length  is  observed  [43].  For  shear 
coaxial  jets,  it  has  been  found  that  the  momentum  flux  ratio  between  the  outer  stream 
and  the  inner  jet  (J  =  (peu2e)  /  (piiif))  is  of  major  importance  for  the  mixing  efficiency 
[22,  69].  It  has  also  been  found  that  combustion  is  more  stable  and  more  efficient  at  high 
values  of  this  quantity  [5].  Different  flow  regimes  may  be  obtained  by  varying  this  ratio 
[51]  and  that  the  inner  jet  core  length  decreases  as  the  momentum  flux  ratio  increases 
[13,  25], 

The  dynamics  of  coaxial  transcritical  jets  is  of  interest  for  the  study  of  combustion  in- 
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stability.  This  can  be  examined  by  considering  the  interaction  of  acoustic  waves  with 
transcritical  flows  and  its  influence  on  mixing  efficiency.  The  effect  of  an  acoustic  modu¬ 
lation  on  a  single  round  jet  at  sub-  and  supercritical  pressure  is  investigated  by  Chehroudi 
and  Talley  [8]  who  find  that  effects  on  the  jet  structure  are  more  profound  in  the  subcriti- 
cal  case.  The  impact  of  the  modulation  is  reduced  as  the  initial  jet  velocity  is  augmented. 
The  impact  of  acoustics  on  a  non-reacting  shear  coaxial  jet  has  been  investigated  exper¬ 
imentally  by  Leyva  and  Rodriguez  [25,  24,  55,  54],  They  have  assessed  the  effects  of 
acoustics  on  the  jet  by  measuring  the  spreading  angle  of  the  inner  and  outer  jets  and  the 
length  of  the  inner  jet  before  its  first  break-up.  It  has  been  found  that  the  injector  geom¬ 
etry  has  an  effect  on  the  susceptibility  of  the  coaxial  flows  to  acoustic  modulation.  For 
the  injector  geometry  studied  here,  acoustics  have  a  measurable  effect  on  the  reduction  of 
the  inner  jet  dark  core  length  for  an  intermediate  value  of  J  at  pressures  varying  from  1.5 
to  5.0  MPa.  Reactive  cases  were  also  experimentally  investigated  (see  [5,  19,  67,  64,  65] 
for  example,  and  [53,  32]  for  acoustically  modulated  cases)  but  will  not  be  detailed  here 
since  they  are  out  of  the  scope  of  the  present  article. 

The  numerical  modeling  of  such  complex  flows  has  been  considered  more  recently  by 
various  groups.  It  is  first  noted  that  in  the  transcritical  range,  thermodynamic  properties 
notably  differ  from  those  of  a  perfect  gas  and  cannot  be  accurately  represented  with  the 
standard  perfect  gas  (PG)  equation  of  state  (EOS).  Thermodynamics  and  “real-gas”  (RG) 
EOS  are  required  [47,  15,  71].  Cubic  equations  of  state  are  generally  used  to  this  purpose 
[37],  but  more  precise  (and  numerically  expensive)  equations  may  also  be  employed  [71]. 
This  issue  and  the  modeling  of  transport  properties  are  extensively  reviewed  in  [1],  [37], 
[31]  and  [2], 

The  models  described  in  the  previous  studies  are  used  in  various  calculations.  Mixing  at 
supercritical  pressure  is  studied  with  direct  numerical  simulation  in  a  mixing  layer  config¬ 
uration,  to  show  the  role  of  density  gradients  on  the  global  layer  stability  and  turbulence 
characteristics  and  the  strong  effect  on  the  molecular  mixing  of  the  reduced  species  dif¬ 
fusion  near  critical  conditions  [2],  This  highlights  the  need  for  approriate  subgrid-scale 


models  for  Large  Eddy  Simulation  (LES)  under  supercritical  pressure  conditions  [63,  68]. 
A  few  large  eddy  simulations  (LES)  have  been  performed  under  transcritical  conditions. 
An  extensive  review  of  non-reacting  LES  in  the  period  ending  in  2006  can  be  found  in 
Zong  et  al.  [76].  A  single  nitrogen  round  jet  was  studied  by  Zong  [74,  73]  and  recently 
by  Schmitt  et  al  [59].  Again,  the  stabilizing  effect  of  the  density  gradient,  and  its  role  on 
turbulent  energy  redistribution  along  the  mixing  layer  was  identified  [74],  Shear  coaxial 
injection  of  oxygen  and  methane  has  also  been  simulated  by  Zong  and  Yang  [75].  This 
has  been  continued  by  considering  a  non- reactive  coaxial  injector  submitted  to  acoustic 
modulations  [26]  where  a  significant  reduction  of  the  jet  intact  core  length  is  observed 
when  acoustic  perturbations  are  present. 

One  also  notes  that  LES  have  been  used  for  reactive  cases  at  supercritical  pressures. 
Shear  coaxial  jet  flames  were  investigated  by  Oefelein  [37,  38],  Oefelein  and  Yang  [36], 
Zong  and  Yang  [77],  Matsuyama  et  al.  [28],  Masquclet  et  al.  [27]  and  Schmitt  et  al. 
[57].  The  problem  is  also  considered  in  the  Reynolds  Average  Navier-Stokes  (RANS) 
framework  (see  for  example  [16],  [11]  or  [49]). 

The  present  investigation  focuses  on  the  mixing  behavior  in  the  flow  originating  from  a 
shear  coaxial  injector  in  the  absence  or  in  the  presence  of  external  acoustic  modulation. 
Three  operating  points  are  considered  corresponding  to  different  momentum  flux  ratios. 
The  coaxial  jet  geometry  corresponds  to  a  thick  inner  jet  post  (the  wall  which  separates 
the  inner  and  outer  streams  in  the  coaxial  configuration  has  a  thickness  of  the  order  of  the 
jet  diameter).  In  this  article  each  case  is  defined  by  the  chamber  pressure,  the  momentum 
flux  ratio,  and  the  coaxial  jet  geometry  which  is  fixed  in  the  present  study.  The  main 
components  of  the  experimental  facility  are  reviewed  in  section  2.  Balance  equations 
and  models  used  in  the  simulations  are  briefly  described  in  section  3.  The  definition  of 
the  computational  domain  and  other  numerical  aspects  are  considered  in  section  4.  The 
influence  of  the  momentum  flux  ratio  is  examined  in  section  5.  Results  obtained  are 
used  to  examine  the  dynamics  of  the  flow  and  identify  the  spectral  content  in  the  various 
regions.  Effects  of  a  transverse  acoustic  modulation  are  investigated  in  section  6. 
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2  Experimental  configuration 


The  experimental  work  described  in  this  article  was  carried  out  at  the  Supercritical 
Cryogenic  Test  Facility  (EC-4)  of  the  Air  Force  Research  Laboratory  at  Edwards  Air 
Force  Base,  California.  The  facility  features  a  test  chamber  that  can  achieve  pressures 
over  5  MPa  and  cryogenic  injection  at  temperatures  below  120  K.  The  photograph  in 
Fig.  1(a)  shows  the  main  test  chamber  and  several  key  components  of  the  facility.  A 
schematic  view  of  the  test  chamber  is  provided  in  Fig.  1(b).  The  facility  was  used  to 
obtain  a  shear  coaxial  jet  flow  where  the  inner  jet  was  cooled  down  to  be  denser  than 
the  outer  jet  under  typical  combustion  chamber  pressures,  resembling  conditions  found 
in  practice  where  the  inner  jet  oxidizer  is  injected  at  a  significantly  lower  temperature 
compared  to  that  of  the  outer  jet  stream.  A  second  equally  important  objective  was  to 
include  the  capability  to  perform  acoustic  forcing.  When  this  was  achieved  the  coaxial  jet 
flow  parameters  were  varied  and  the  coaxial  jet  behavior  was  observed  with  and  without 
acoustic  excitation. 


Heat 

exchanger 

Coaxial 

/x/injector  Acoustic 
i*  l  driver  unit 

(piezo-siren) 


Chamber  pressurization 
Chamber  vent 


Figure  1:  (a)  Photograph  of  the  main  test  chamber  and  key  instrumentation  of  EC4, 
the  supercritical  cryogenic  facility  at  AFRL,  Edwards  AFB.  (b)  Flow  diagram  of  the  test 
chamber  and  its  vicinity. 


To  run  a  test,  a  laboratory-wide  high  pressure  gaseous  nitrogen  supply  was  fed  into  the 
inner  and  outer  flow  lines  of  the  coaxial  jet  and  the  flow  rates  adjusted  to  correspond  to 
those  of  the  desired  test  conditions.  The  mass  flow  rates  through  the  inner  and  outer  jets 
were  measured  with  Porter  mass  flow  meters  (models  122  and  123-DKASVDAA).  The 
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same  gaseous  nitrogen  supply  used  for  the  coaxial  jet  was  also  employed  to  increase  the 
pressure  in  the  test  chamber  and  achieve  the  required  level,  which  was  measured  with  a 
Stellar  1500  transducer.  The  jet  temperatures  were  controlled  by  adjusting  the  flow  rates 
of  liquid  nitrogen  through  the  heat  exchangers  placed  on  the  feed  lines.  Densities  of  the 
two  jets  were  deduced  from  the  measured  temperatures  and  pressure  using  the  National 
Institute  of  Standards  and  Technology’s  thermophysical  properties  of  fluid  systems  online 
database  [23].  These  values  in  conjunction  with  the  mass  flow  rate  data  provided  the 
injection  velocities  and  the  corresponding  Reynolds  numbers,  velocity  and  momentum 
flux  ratios.  Several  thermocouples  were  used  across  the  heat  exchanger  and  in  other 
locations  to  keep  track  of  the  conditions  of  the  flow  in  order  to  maintain  the  required  flow 
properties  for  each  test  case.  Appendix  A  provides  a  summary  of  the  operating  conditions 
achieved  for  all  cases  reported  in  the  present  article.  Two  images  of  the  shear  coaxial 
injector  can  be  seen  in  Fig.  2.  The  first,  Fig.  2(a),  gives  a  full  view  before  assembly  and 
the  second,  Fig.  2(b)  shows  a  plane  view  of  the  end  section.  The  inner  jet  diameter  is 
di  =  0.51  mm.  The  outer  jet  has  an  inner  diameter  of  1.59  mm  and  outer  diameter  of 
2.42  mm.  The  injector  outer  diameter  was  3.18  mm.  The  length  to  inner  diameter  ratio 
was  100  for  the  inner  jet  and  34  for  the  outer  jet  based  on  the  hydraulic  diameter.  The 
coaxial  injector  was  installed  so  that  the  inner  and  outer  jets  were  nominally  concentric. 
The  inner  jet  exit  plane  was  recessed  by  0.3  mm  from  the  outer  jet  exit  plane.  This  recess 
length  was  chosen  to  mimic  realistic  coaxial  jet  configurations  used  in  practice. 

The  coaxial  injector  was  placed  in  an  inner  chamber  within  the  main  test  chamber.  This 
unit  was  protruding  into  the  chamber  by  5  mm  through  a  17  mm  in  diameter  central 
orifice  in  the  top  wall.  This  inner  chamber  was  created  to  augment  the  level  of  acoustic 
modulation  acting  on  the  coaxial  flow.  The  inner  chamber  had  a  height  of  6.6  cm,  a 
width  of  7.6  cm  and  a  depth  of  1.3  cm.  The  flow  was  exhausted  through  another  orifice 
at  the  center  of  the  bottom  wall. 

Injection  temperatures  were  measured  with  unshielded  traveling  type  E  thermocouples 
with  a  bead  diameter  of  0.1  mm.  The  accuracy  of  the  thermocouples  used  in  the  study  was 
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* 


Figure  2:  (a)  Photograph  of  the  shear  coaxial  injector  used  in  this  study,  (b)  Exit  plane 
view  of  the  injector. 


checked  with  an  RTD  and  found  to  be  within  1  K.  A  miniature  pressure  transducer  (Kulite 
Semiconductor  Products,  models  CCQ-062-1000A  and  CCQ-093-750A)  was  placed  next 
to  the  thermocouple  using  a  small  metallic  post  for  support,  and  recorded  pressure  fluctu¬ 
ations  at  a  sampling  frequency  of  20  kHz.  These  pressure  transducers  having  an  absolute 
pressure  range  of  either  6.9  MPa  or  5.2  MPa,  respectively  were  used  to  used  to  scan  the 
waveheld  induced  by  the  acoustic  driver  units  prior  to  flow  experiments. 

The  coaxial  jet  structure  was  observed  with  backlighting  (shadowgraph)  imaging.  The 
flow  was  illuminated  with  a  light  source  (Newport  variable  power  arc  lamp  set  at  160W  or 
300W)  and  the  transmitted  light  was  recorded  with  a  Phantom  v7.1  high-speed  camera 
capable  of  recording  up  to  160,000  frames  per  second  at  a  resolution  of  32  pixels  by  32 
pixels.  The  camera  equipped  with  AF  Nikkor  35-105  mm  lens  and  a  Nikon  No.  1  Close- 
Up  lenscan  be  seen  at  the  bottom  center  in  Fig.  1(a).  The  image  resolution  in  the  present 
experiments  varied  from  128  pixels  x  224  pixels  to  196  pixels  x  400  pixels,  depending 
on  chamber  pressure  and  outer-to-inner  momentum  flux  ratio.  The  images  were  in  a 
grayscale  scale  with  each  pixel  taking  a  value  between  0  and  255  for  different  intensities 
of  gray  ranging  from  white  to  black.  Each  pixel  represents  an  approximate  area  of  0.08 
mm  x  0.08  mm.  Depending  on  the  resolution  chosen,  the  framing  rate  was  20,  25  or  41 
kHz.  The  exposure  time  was  generally  1-2  /i s  and  the  number  of  images  saved  per  run 
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was  1000. 


The  images  recorded  were  then  processed  using  a  Matlab  subroutine  based  on  the  Otsu 
technique  [44] .  The  subroutine  takes  a  grayscale  image  and  finds  a  pixel  threshold  value. 
All  the  pixel  values  above  this  threshold  are  converted  to  black  and  the  remaining  pixels 
to  white.  Then  the  subroutine  measures  the  length  of  the  black  region  extending  from  the 
exit  of  the  coaxial  jet.  This  is  referred  to  as  the  “dark  core  length”.  A  “mean  dark  core 
length”  lac  is  then  obtained  by  averaging  over  the  set  of  images  recorded  by  the  camera. 
Acoustic  modulation  of  this  setup  was  obtained  with  two  piezo-sirens  placed  at  each  end 
of  the  chamber  at  a  distance  of  34  cm  from  the  jet  axis,  custom-designed  for  the  Air 
Force  Research  Laboratory  by  Hersh  Acoustical  Engineering,  Inc  (the  piezo-sirens  can  be 
seen  at  both  sides  of  the  chamber  in  Fig.  1(a)).  A  Fluke  100  MS/s  arbitrary  waveform 
generator  (model  292)  was  used  to  produce  two  sinusoidal  waves  with  the  same  frequency 
but  with  a  prescribed  phase  between  them.  The  signals  were  then  sent  to  two  amplifiers 
(Krohn-Hite  model  7500  and  a  Trek  model  PZD2000A),  one  for  each  piezo-siren.  The 
amplified  signals  were  in  the  200  to  540  V  range.  To  accommodate  for  the  rectangular 
chamber,  a  waveguide  with  a  catenary  contour  -  see  upper  left  corner  of  Fig.  1(a)  -  was 
used  to  transmit  the  waves  from  the  circular  cross-section  of  each  siren  to  a  rectangular 
cross-section.  The  maximum  root-mean-square  acoustic  pressure  fluctuations  generated 
by  the  two  piezo-ceramic  acoustic  sources  in  the  inner  chamber,  measured  using  the  three 
differential  Kulite  pressure  transducers  (model  XCQ-093-25D)  flush  mounted  on  a  lateral 
chamber  wall,  varied  from  8.0  to  22  kPa  and  were  produced  in  the  2.9  to  3.1  kHz  range. 
These  operating  conditions  were  obtained  by  manually  varying  the  frequency  on  the  signal 
generator  and  finding  the  highest  possible  pressure  amplitude. 

The  phase  difference  between  the  signals  sent  to  the  piezo-siren  elements  was  varied  to 
expose  the  coaxial  jet  to  different  effective  positions  relative  to  the  pressure  node  or 
antinode  of  the  acoustic  field.  When  the  two  piezo-siren  elements  produce  waves  with  a 
zero  degree  phase  angle  between  them,  the  motion  of  the  piezo-siren  radiating  membranes 
is  synchronized  and  radiates  in  opposite  directions.  This  produces  conditions  of  high 
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pressure  perturbations  and  low  velocity  fluctuations  at  the  center  of  the  chamber,  where 
the  coaxial  jet  is  located,  which  correspond  to  a  pressure  antinode  or  velocity  node.  In 
contrast,  when  the  two  drivers  present  a  180-degree  phase  difference,  then  the  radiating 
membranes  move  in  the  same  direction,  generating  high  velocity  fluctuations  with  very 
small  variations  in  pressure  at  the  jet  location,  producing  a  pressure  node  or  velocity 
antinode. 


3  Balance  equations  and  models 

At  this  stage  it  is  worth  briefly  reviewing  the  spatially  filtered  balance  equations  and 
introducing  the  real  gas  models  and  subgrid  scale  closures  used  in  the  large  eddy  simula¬ 
tions. 


3.1  Governing  equations 


A  single  species  flow  configuration  is  investigated  in  the  present  study.  In  this  situation, 
the  vector  of  conservative  variables  for  a  compressible  flow  is  w  =  ( p ,  pu,  pE)T ,  where  p  is 
the  density,  u  the  velocity  vector  and  E  the  total  energy  ( E  is  the  sum  of  internal  energy 
es  and  kinetic  energy  e*,  =  1/2^-uf).  The  present  large-eddy  simulations  are  carried 
out  by  integrating  the  mass  weighted  spatially  filtered  Navier-Stokes  equations.  In  the 
following  balance  equations  (for  a  given  variable  0,  0  designates  the  spatially  filtered 
variable,  while  0  is  the  mass-weighted  (Favre)  spatially  filtered  variable)  : 


dp 

dt 


+ 


dpuj 

dxj 


dpUi 

dt  + 


dpujiii 

dxj 


dpE  dpEuj 
dt  dxj 


dpuj 

dxj 


dp  dr ij_  _  drtj 

dxi  dxj  dxj 

dq ]  _  dqj_  +  dupr^ 
dxj  dxj  dxj 


(1) 

(2) 
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In  Eqs.  1-3,  p  is  the  pressure,  x  the  spatial  coordinates  vector  and  t  is  the  time.  The 
laminar  viscous  stress  tensor  r  and  heat  flux  vector  q  are  expressed  as  linear  functions 
of  the  strain  rates  and  temperature  gradient  respectively.  Laminar  viscosity  and  heat 
conductivity  coefficients  are  determined  with  the  Chung  et  al.  method  [12].  For  the 
single  species  case  investigated  in  the  present  article,  no  additional  diffusion  terms  are 
needed  since  the  Soret  and  Dufour  effects  naturally  vanish.  The  modeling  of  the  subgrid- 
scale  stress  tensor  rh  and  heat  flux  vector  ql  requires  some  discussion.  It  is  known  from 
recent  studies  that  subgrid  closure  rules  for  real  gases  differ  from  those  used  for  perfect 
gas  flows.  In  principle  one  should  take  into  account  additional  terms  arising  from  the 
nonlinearity  of  the  state  equation,  and  from  the  large  variations  in  transport  properties 
as  explained  in  some  recent  studies  by  Bellan  et  al.  [2],  Scllc  et  al.  [63]  and  Taskinoglu 
and  Bellan  [68].  However  the  modeling  issues  are  not  completely  settled.  It  is  also  known 
from  other  calculations  that  the  low  pressure  subgrid  scale  models  yield  suitable  results 
for  transcritical  and  supercritical  round  jets  [59,  62],  It  was  decided  for  these  reasons  to 
keep  these  low  pressure  closures  for  the  subgrid  scale  models  as  described  in  more  detail 
in  Sec.  3.3. 


3.2  Real  gas  thermodynamics  and  equation  of  state 

Real  gas  thermodynamics  for  high  pressure  flow  calculations  is  envisaged  in  many  recent 
studies  where  accurate  descriptions  of  the  fluid  state  are  derived  from  various  types  of 
equations  [1,  71].  The  present  simulations  rely  on  the  Peng- Robinson  equation  of  state 
[45]  (Eq.  4)  which  is  used  in  the  analysis  to  model  the  departure  from  an  ideal-gas 
behavior.  This  EOS  is  less  accurate  than  the  modified  version  but  it  is  also  simpler  and 
requires  a  smaller  amount  of  calculations.  The  fluid  state  is  expressed  in  the  form  : 


prT  p2a(T) 

1  —  pb  1  +  2  pb  —  p2b2 


(4) 
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where  T  is  the  temperature,  r  =  R/W  where  R  is  the  perfect-gas  constant  and  W  the 
molar  mass.  The  Peng- Robinson  coefficients  a(T)  and  b  for  a  single-species  fluid  are  [47]: 


a  =  0.457236-^—^-  (l  +  c(l  —  a/TVjV  (5) 

Pc  \  / 

tT 

b  =  0.077796—  (6) 

Pc 


with 


c  =  0.37464  +  1.54226c;-  0.26992c;2 


(7) 


where  Tc  =  126.2  K  and  pc  =  3.396  MPa  are  the  critical  temperature  and  pressure  of 
nitrogen,  u  =  0.0372  represents  its  acentric  factor  and  Tr  =  T/Tc  is  the  reduced  temper¬ 
ature.  This  equation  offers  a  good  trade-off  between  computational  cost  and  precision. 
The  pressure  dependence  of  thermodynamic  functions  and  coefficients  is  subsequently  de¬ 
rived  from  the  EOS  [33,  38,  71].  Departure  energy  Aes  and  specific  heat  Ac„  at  constant 
volume  are  given  by  [20,  47]: 


Ae.s 


A  cv 


(8) 

(9) 


where  Aes  =  es(p,  T )  —  e°(T)  and  A cv  =  cv(p,  T )  —  c°(T).  The  enthalpy  and  specific  heat 
at  constant  pressure  are  then  deduced  from  : 
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The  thermal  expansion  coefficient  a  and  the  isothermal  compressibility  coefficient  (3  ap¬ 
pearing  in  the  previous  expressions  are  given  by: 


a  =  — 


1 

P 


p 


P 


1  /  dp\ 

P  \dp) T 


(11) 


The  low-pressure  references  (i.e.  e°  and  c°)  are  obtained  from  the  JANAF  Thermody¬ 
namical  Tables  [7],  and  only  depend  on  temperature. 


3.3  Subgrid-scale  models  and  assumptions 

As  indicated  previously,  the  present  simulations  rely  on  a  standard  (low  pressure)  subgrid 
scale  closure,  the  WALE  model,  in  which  the  stress  tensor  rl  is  conveniently  expressed 
in  terms  of  the  spatially  resolved  strain  rates  [35]  : 


v 


2i/ 


(12) 


where 


S  =  (C„A)' 


(44)3/2 


(SySy)^  +  (44)5/4 


sij 
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(13) 

(14) 


In  the  previous  expression  the  WALE  model  constant  Cw  =  0.4929  and  gij  is  the  irrota- 
tional  part  of  the  stress  tensor  : 


9ij 


duj 

dxj 


and  fpj  (Jikgkj 


(15) 


The  WALE  model  distinguishes  velocity  gradients  associated  with  the  rotational  and 
pure  shear  components  of  the  velocity  gradients  in  order  to  distinguish  regions  where 
turbulence  is  transitional  and  regions  where  it  is  is  fully  developed.  It  is  well  suited  to 
the  treatment  of  shear  flows  of  the  type  considered  in  this  article. 
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The  subgrid  scale  heat  flux  is  modeled  using  a  standard  gradient  transport  assumption  : 


ft*  -A^  (16) 

The  subgrid  scale  conductivity  is  cast  in  the  form  Xt  =  cppt/¥xt  where  fit  =  pvt  and  Prt 
respectively  designate  the  subgrid  dynamic  viscosity  and  turbulent  Prandtl  number.  As 
in  low  pressure  simulations  a  constant  value  Pr*  =  0.7  is  used  in  the  present  calculations. 

4  Numerical  aspects 

The  mass  weighted  spatially  filtered  compressible  Navier-Stokes  equations  in  combination 
with  the  subgrid  scale  models  described  previously  are  integrated  in  the  AVBP  flow 
solver  [60,  34],  This  code  has  already  been  used  to  calculate  low  and  moderate  pressure 
reactive  flows  of  multicomponent  mixtures  of  perfect  gases.  The  spatial  discretization  on 
unstructured  or  hybrid  meshes  facilitates  applications  in  which  the  geometry  is  complex, 
a  feature  which  is  quite  useful  if  one  wishes  to  analyze  practical  configurations.  The  LES 
of  turbulent  flows  requires  a  low-dissipation  algorithm  [21,  18]  which  is  often  obtained 
with  centered  discretization  schemes.  The  present  integration  method  relies  on  a  Taylor- 
Galerkin  weighted  residual  central  distribution  scheme,  called  TTG4A.  This  scheme  is 
third-order  in  space  and  fourth-order  in  time  [50]. 

The  real  gas  equation  of  state,  thermodynamics  and  transport  coefficients  have  been  im¬ 
plemented  in  AVBP.  The  convective  fluxes  Jacobian  matrices,  used  by  the  scheme,  are 
expressed  in  terms  of  real  gas  thermodynamics  to  preserve  the  overall  consistency  of  the 
code.  A  fully  consistent  treatment  of  associated  boundary  conditions  is  based  on  the 
characteristic  wave  decomposition  method  NSCBC  [46,  34]  and  includes  specific  expres¬ 
sions  provided  for  real  gas  thermodynamics  in [39].  The  highly  nonlinear  thermodynamics 
of  the  transcritical  fluid  stream  induces  large  density  gradients  between  the  dense  tran- 
scritical  fluid  and  the  surrounding  gaseous  stream  which  require  a  specific  stabilization 
procedure.  This  relies  on  artificial  viscosity  and  is  used  when  unresolved  gradients  are 
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detected  following  a  method  described  in  [59] .  The  real  gas  AVBP  flow  solver  has  already 
been  used  to  simulate  transcritical  cryogenic  flames  of  various  types  [48,  58,  57]. 


4.1  Computational  domain  and  mesh 


The  geometry  of  the  computational  domain  essentially  matches  the  experimental  config¬ 
uration.  A  three-dimensional  view  of  the  external  boundaries  is  given  in  Fig.  3(a).  The 
longitudinal  section  of  the  chamber  is  rectangular  (59.4x76.0  mm)  and  its  width  is  12.7 
mm.  The  fluid  streams  are  injected  by  a  coaxial  element  which  protrudes  by  5  mm  into 
the  chamber  (Fig.  3b).  This  injector  is  located  in  the  center  of  a  17  nun-diameter  hole  in 
the  inlet  wall.  No  co-flow  is  directly  injected  through  this  orifice  during  the  experiment, 
but  nitrogen  gas  may  traverse  this  aperture  as  a  result  of  entrainment  by  the  coaxial 
streams  giving  rise  to  a  finite  flow  velocity  in  this  region.  The  entrained  fluid  originates 
from  a  large  reservoir  defined  on  the  upstream  side  of  the  central  hole.  Details  about 
the  injector  geometry  are  given  in  Fig.  4(a).  The  inner  injector  diameter,  inner  duct 
thickness,  annular  duct  thickness  and  outer  duct  diameter  are  respectively  d*  =  0.51  mm, 
k  =  0.54  mm,  he  =  0.415  mm  and  dQ  =  2.42  mm.  The  mesh  comprises  2  100  000  nodes 


59.4  mm 


Figure  3:  (a)  Three-dimensional  visualization  of  the  reservoir,  (b)  Longitudinal  cut  of 
the  mesh. 


corresponding  to  10  000  000  tetrahedra.  It  is  highly  refined  near  the  injector  (Fig.  4(b)), 
with  a  constant  characteristic  size  of  0.032  mm  on  a  distance  of  10  inner  jet  diameters. 
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The  mesh  is  then  slowly  coarsened  towards  the  domain  outlet  (Fig.  3(b)). 


Figure  4:  (a)  Closer  view  of  the  injector,  (b)  Longitudinal  cut  of  the  mesh  (35  inner 
injector  diameters). 


4.2  Injection  characteristics  and  boundary  conditions 

Three  test  conditions  are  considered  in  this  analysis  designated  as  “N2”,“N6”  and  “N8”. 
Case  “N6”  serves  as  a  reference  simulation.  The  two  other  cases,  designated  as  “N2”  and 
“N8”,  are  also  calculated  in  order  to  analyze  effects  of  the  momentum  flux  ratio  on  the 
flow  behavior  (Sec.  5).  The  three  cases  differ  by  their  injection  temperatures  and  veloc¬ 
ities,  yielding  momentum  flux  ratios  J  of  1,  2.6  and  9.1  for  cases  “N2”,  “N6”  and  “N8” 
respectively.  The  corresponding  density  ratios  S  are  5.84,  5.74  and  3.22  respectively. 
Injection  conditions  are  gathered  in  Tab.  12  and  Tab.  13  which  respectively  correspond 
to  the  outer  and  inner  streams.  Injection  conditions  are  placed  on  thermodynamic  plots 
on  Fig.  5.  The  density  and  volumetric  internal  energy  notably  depart  from  a  perfect  gas 
behavior  when  temperature  is  close  and  below  the  critical  value.  The  influence  of  acoustic 
perturbations  will  only  be  investigated  in  a  single  case  corresponding  to“N6”  (Sec.  6). 

In  all  cases,  pressure  is  maintained  at  the  outlet  boundary  using  non-reflecting  charac¬ 
teristic  conditions  at  an  adapted  level  of  3.56  MPa,  and  the  chamber  and  reservoir  are 
initially  filled  with  nitrogen  at  the  same  pressure  and  a  temperature  213  K  (the  corre¬ 
sponding  density  is  59.6  kg  m~3).  All  solid  boundaries  are  treated  as  adiabatic  slip  walls. 
Velocity  perturbations  (with  an  amplitude  of  3  %  of  the  mean  flow)  are  added  at  injec- 
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tion  in  order  to  represent  effects  of  turbulent  fluctuations  using  a  procedure  described  in 


[6,  66,  17]. 


Figure  5:  Thermodynamic  state  calculated  with  the  PR  EOS  ( — )  and  the  perfect  gas 
EOS  (■  ■  •).  (a)  Density  in  terms  of  temperature,  (b)  Volumetric  internal  energy  (he. 
pes)  in  terms  of  temperature.  The  pressure  is  constant  and  equal  to  3.56  MPa.  Symbols: 
o  inner  jet  conditions;  x  outer  stream  conditions;  •  chamber  conditions. 


Case 

J  =  (peU2e)/(piU2i) 

M  =  (peue)/(piUi) 

S  pi/  Pe 

U  ue/ Ui 

N2 

1.05 

0.42 

5.86 

2.48 

N6 

3.05 

0.72 

5.65 

4.15 

N8 

9.3 

1.70 

3.22 

5.48 

Table  1:  Injection  characteristics  of  the  simulated  cases.  J  is  the  momentum  flux  ratio 
between  the  outer  and  the  inner  jets,  M  the  mass  flux  ratio,  S  the  density  ratio  and  U  is 
the  velocity  ratio,  p*,  pe,  Ui  and  ue  are  the  inner  and  outer  jet  densities  and  velocities. 


4.3  Simulation  procedure 

Simulations  are  carried  out  on  a  coarse  mesh,  with  an  initial  solution  corresponding  to 
a  uniform  fluid  at  rest.  The  flow  is  then  qualitatively  established  after  three  convective 
periods.  Results  are  then  interpolated  on  the  final  mesh.  Data  averaging  can  be  started 
after  one  transient  period.  The  simulated  physical  time  is  given  in  Tab.  2  together  with 
the  number  of  simulated  convective  periods.  Finally,  the  CPU  time  (equal  to  the  product 
iT-core  x  Uhours)  needed  to  obtain  the  average  solution  (which  does  not  include  the  CPU 
time  needed  for  the  transient  flow  and  the  simulation  on  the  coarse  mesh)  is  also  given 
in  this  table. 
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Case 

A  ta 

ms 

Ata/rUi 

Ata/rUe 

CPU  time  [h] 

N2 

14.3 

1.9 

4.8 

35  000 

N6 

22.7 

4.4 

18.2 

45  000 

N8 

10.4 

4.0 

21.2 

20  000 

Tabic  2:  Averaging  time  for  cases  “N2”,  “N6”  and  “N8”.  A ta  is  the  physical  averaging 
time.  rUi  represents  one  convective  time  over  35  dt  for  the  inner  jet  and  rUe  represents 
one  convective  time  over  35  d,  for  the  outer  jet.  rUi  =  (35 di)/ui  ;  rUe  =  (35 di)/ue. 

5  Effect  of  momentum  flux  ratio 


It  is  first  natural  to  consider  the  influence  of  the  momentum  flux  ratio  on  the  system 
behavior.  Flow  patterns  calculated  numerically  and  obtained  experimentally  are  first 
compared  in  Sec.  5.1.  This  is  followed  by  a  description  of  the  mean  flow  in  Sec.  5.2  and 
by  an  analysis  of  flow  dynamics  in  Sec.  5.3. 


5.1  Instantaneous  flow  patterns 

Instantaneous  backlighting  images  obtained  experimentally  are  shown  in  Fig.  6(a),(c),(e). 
In  the  dark  core  light  is  deviated  by  changes  in  the  refraction  index.  This  in  turn  delin¬ 
eates  regions  where  the  density  takes  large  values  (the  deviation  can  be  roughly  linked 
to  the  second  derivative  of  the  density  integrated  on  the  line  of  sight).  It  is  generally 
admitted  that  the  dark  region  represents  the  position  of  the  dense  core,  and  qualitatively 
highlights  its  external  boundaries  giving  an  insight  on  the  main  structures  present  in  the 
flow  [43].  As  the  momentum  flux  ratio  is  increased,  the  inner  jet  axial  length  is  reduced. 
Case  “N2”  exhibits  large  scale  structures  over  about  10  inner  injector  diameters  down¬ 
stream  of  the  injector.  These  structures  are  not  identifiable  in  cases  “N6”  and  “N8”. 
The  latter  is  characterized  by  a  sudden  termination  of  the  inner  jet  at  a  few  inner  jet 
diameters  from  the  exit  plane.  In  contrast  with  cases  “N2”  and  “N6”,  the  inner  stream 
appears  to  be  quickly  mixed  with  the  surrounding  fluid,  and  the  central  jet  can  no  longer 
be  clearly  identified  downstream  of  this  abrupt  termination. 

Longitudinal  slices  of  instantaneous  density  distributions  extracted  from  simulations  are 
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shown  in  Fig.  6(b) ,  (d)  ,(f )  which  respectively  correspond  to  “N2”,  “N6”  and  “N8”  cases. 
Comparisons  between  these  figures  and  the  experimental  backlighting  images  are  ad¬ 
mittedly  qualitative.  Nevertheless,  general  experimental  features  corresponding  to  these 
different  cases  are  replicated.  One  finds  that  case  “N2”  features  large  scale  structures 
which  are  not  present  in  cases  “N6”  and  “N8”.  The  abrupt  termination  noticed  in  case 
“N8”  is  retrieved  numerically.  This  is  also  made  more  apparent  in  the  density  iso-surfaces 
corresponding  to  p  —  0.1  (p*  —  pe)  +  pe )  plotted  in  Fig.  7.  The  large  scale  structures  ob¬ 
served  experimentally  characterize  case  “N2”,  and  disappear  as  the  momentum  flux  ratio 
is  increased  (cases  “N6”  and  “N8”). 


5.2  Mean  flow  description 

It  is  first  interesting  to  examine  effects  of  the  momentum  flux  ratio  on  the  mean  “dark 
core”  length.  Since  the  density  field  cannot  be  directly  compared  to  the  experimental 
light  distribution  images,  it  is  natural  to  define  a  set  of  density  length  scales  and  compare 
them  with  the  dark  core  length  deduced  from  the  data.  The  jet  core  characteristic  lengths 
are  defined  by  considering  the  axial  density  distribution  : 


/0.99 

V 


<0.5 

V 


-  x0 


lO.l 

V 


-  x0 


(17) 


where  Xo  is  the  axial  position  of  the  inner  injector  exit,  x°p",  x°p 5  and  x J)'1  are  positions  on 
the  jet  axis  where  density  is  respectively  equal  to  99  %  of  its  initial  value,  to  0.5 (pi—pe)+pe 
or  to  0.1  (pi  —  pe)  +  pe.  These  three  lengths  are  shown  in  Fig.  8(a)  together  with  mean 
density  profiles  on  the  centerline. 

It  is  also  useful  to  define  at  this  point  another  length  scale  based  on  the  axial  mean 
velocity  profile.  The  central  jets  all  feature  an  initial  decrease  in  the  axial  velocity  on  the 
centerline  as  can  be  seen  in  Fig.  8(b).  The  minimum  value  is  reached  at  xu  which  is  used 
to  define  a  new  length  scale  lu  =  xu  —  x0. 

Characteristic  lengths  are  gathered  in  Tab.  3  for  comparison  with  experimental  mea- 
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□ 

(f) 


Figure  6:  Comparison  between  experimental  visualizations  (backlighting)  and  typical  in¬ 
stantaneous  density  distributions  obtained  from  simulations  (white:  60  kg  m”3  ;  black: 
inner  jet  injection  density  (Tab.  13)  ;  logarithmic  scale).  Case  “N2”  (J=1.05)  :  (a) 
Experimental  backlighting  image,  (b)  Calculated  density.  Case  “N6”  (J=3.05)  :  (c) 
Experimental  backlighting  image,  (d)  Calculated  density.  Case  “N8”  (J=9.3)  :  (e)  Ex¬ 
perimental  backlighting  image,  (f)  Calculated  density. 
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4 


□ 


(b) 


□ 


□ 

(c) 

Figure  7:  LES  results.  Density  iso-surface  corresponding  to  p  —  0.1  (pi  —  pe )  +  pe  colored 
by  the  axial  velocity,  (a)  case  “N2”  (J=1.05)  umin=- 2  m  s^2,  umax= 6.5  m  s-2,  (b)  case 
“N6”  (J =3.05)  umin=-4:  m  s-2,  umax= 16  m  s'2,  (c)  case  “N8”  (J =9.3)  umin=- 5  m  s-2, 
umax= 40  m  s'2.  The  plot  corresponds  to  20  inner  injector  diameters. 
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surements  of  the  “dark  core”  length.  The  latter  values  are  deduced  from  backlighting 
images  by  measuring  the  “dark  region”  size  from  average  images.  Again,  comparison 
between  experiment  and  simulation  is  limited.  As  already  noted  in  [43] ,  lengths  deduced 
from  experimental  backlighting  data  are  generally  larger  than  those  deduced  from  density 
measurements.  It  is  however  clear  that  there  are  similarities  between  trends  in  simulations 


and  experiments.  In  agreement  with  experiments,  numerical  results  indicate  a  reduction 
in  the  density  length  scales  as  the  momentum  flux  ratio  is  increased.  The  velocity  length 
scale  is  reduced  and  the  amplitude  of  the  velocity  defect  is  enhanced  as  the  momentum 
flux  ratio  is  increased.  The  flow  may  also  be  characterized  by  the  distributions  of  mean 


Figure  8:  Centerline  profiles  deduced  by  averaging  LES  results,  (a)  Density,  (b)  Axial 
velocity  normalized  by  the  outer  injection  velocity  u*  =  u/ue.  —  :  Case  “N2”,  —  •  —  : 
Case  “N6”  and - :  Case  “N8”. 


Case 

IflAi 

lu/  di 

Idc/  di 

N2 

2.57 

4.40 

8.66 

3.8 

15.5 

N6 

1.81 

3.37 

5.88 

2.95 

6.0 

N8 

0.85 

1.17 

1.67 

1.96 

2.5 

Table  3:  Characteristic  lengths  deduced  from  simulations  based  on  axial  density  and 
velocity  profiles  and  from  experimental  measurement  of  the  dark  core  size. 


axial  velocity  plotted  in  Fig.  9.  The  outer  flow  contracts  towards  the  central  stream 
and  pinches  the  central  jet  to  a  point  where  this  stream  disappears.  In  the  three  cases, 
a  small  back-flow  is  generated  behind  the  injector  lips.  As  the  momentum  flux  ratio  is 
increased  beyond  a  certain  critical  value  (case  “N8”),  the  central  jet  is  surrounded  by  a 
negative  axial  velocity  region  (Fig.  9c). 
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(a) 


(b) 


(c) 

Figure  9:  Distribution  of  axial  velocity  (white:  minimum  ;  black:  maximum).  This  map 
corresponds  to  35  inner  injector  diameters.  —  iso-contour  of  zero  axial  velocity,  (a)  case 
“N2”  (J—1.05),  (b)  case  “N6”  ( J=3.05 ),  (c)  case  “N8”  (J=9.3). 
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Case  “N8”  exhibits  a  large  negative  value,  as  observed  previously  in  Fig.  9(c).  Further 
downstream,  the  flows  behave  in  a  self-similar  fashion  in  terms  of  axial  velocity  normalized 
by  the  outer  jet  injection  velocity.  Beyond  this  section,  the  impact  of  the  inner  jet  on 
the  flow  disappears  as  this  jet  is  mixed  with  the  outer  stream.  Since  the  mass  flow  rate 
of  the  outer  jet  is  much  higher  than  that  of  the  inner  jet,  the  outer  jet  dominates  the 
flow.  Interestingly,  the  three  profiles  merge  at  x  >  15 di,  indicating  a  self-similar  behavior 
when  profiles  are  plotted  in  terms  of  normalized  variables. 

Radial  profiles  of  normalized  density  and  axial  velocity  are  plotted  in  Fig.  10  for  case 
“N6”  (case  “N2”  is  not  shown  here  because  its  behavior  is  quite  similar).  The  first 
profiles  lie  in  the  potential  core  region  (x  <  x°p"),  and  feature  constant  levels  of  density 
and  velocity  in  that  region.  Negative  values  of  axial  velocity  indicate  that  back-flow  is 
generated  behind  the  lips.  The  density  profile  behind  the  inner  lip  shows  a  rapid  variation 
towards  the  outer  stream.  This  is  attributed  to  an  efficient  mixing  of  inner  and  outer 
streams  associated  with  the  formation  of  vortical  structure  behind  the  inner  lip  and  the 
generation  of  a  back  flow.  In  this  region,  mass  is  extracted  from  the  central  jet  by  vortical 
structures  and  recirculates  upstream.  Dynamics  in  this  region  are  discussed  in  Sec.  5.3. 
The  profile  at  3d*  corresponds  to  x  =  xu.  The  density  stratification  is  still  apparent.  At 
x  =  15 di,  velocity  and  density  are  maximum  on  the  centerline,  and  a  jet-like  profile  is 
retrieved.  Further  downstream,  the  flow  behaves  like  a  single  round  jet  and  develops  in  a 
self-similar  fashion.  The  half  width  at  half  maximum  (HWHM)  of  density  Lp  is  plotted 
in  Fig.  11(a).  For  x  >  20 di,  this  quantity  evolves  quasi-linearly  Lp  =  Ep(x  —  xv),  where 
E h  is  the  spreading  rate  of  the  jet  and  xv  is  a  virtual  origin.  Linear  regressions  for 
20 di  <x<  35dj  give  spreading  rates  of  0.12,  0.11  and  0.11  for  cases  “N2”,  “N6”  and  “N8”, 
respectively.  These  values  are  in  agreement  with  those  given  for  example  by  Chen  and 
Rocli  [10]  for  a  single  jet  and  with  results  obtained  from  previous  high  pressure  simulations 
of  single  species  round  jets  [59,  76].  In  this  region  of  the  flow,  the  density  profile  takes  a 
gaussian  shape  (Fig.  lib).  A  linear  evolution  of  the  density  decay  is  also  observed,  but 
not  shown  here,  confirming  that  the  flow  is  self-similar  in  the  downstream  region.  Since 


Figure  10:  LES  results,  Case  “N6”.  (a)  Radial  profiles  of  normalized  density  p*  = 

(p  —  Poo)/{pc  —  Poo )  (b)  Radial  profiles  of  axial  velocity. 


Figure  11:  LES  results  (a)  Half  Width  Half  Max  of  density  —  Case  “N2”,  —  •  —  Case 
“N6”  and - Case  “N8”.  (b)  Radial  profiles  of  normalized  density  p'  =  (p  —  Poo)/(pc  — 

Poo)- 
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the  flow  develops  in  a  confined  domain,  one  may  wonder  if  it  is  influenced  by  recirculation 
in  the  chamber.  Examination  of  velocity  distributions  in  the  domain  indicates  that  the 
outer  recirculation  velocities  are  low  and  that  the  reverse  flow  in  the  outer  region  has  a 
limited  influence  on  the  coaxial  jet  structure. 


5.3  Near  injector  flow  dynamics 

Large  eddy  simulations  allow  detailed  analysis  of  the  flow  dynamics.  This  can  be  accom¬ 
plished  by  estimating  the  spectral  content  of  fluctuations  in  various  regions  of  the  flow  to 
identify  the  preferred  modes  of  the  system.  To  this  end,  it  is  convenient  to  place  a  set  of 
sensors  in  the  flow  and  examine  the  corresponding  velocity  and  density  time  series.  This 
is  illustrated  in  Fig.  12.  Since  the  general  flow  structures  of  cases  “N2”  and  “N6”  differ 
from  that  corresponding  to  “N8”,  it  is  natural  to  study  these  two  cases  separately  (Sec. 
5.3.1).  Case  “N8”  is  then  considered  in  Sec.  5.3.2. 


(a) 


0  1  2  3  4  5  £  \ 
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Figure  12:  The  three  sensor  locations.  The  numerical  probes  are  used  to  measure  char¬ 
acteristic  frequencies  in  the  coaxial  flow,  (a)  Case  “N6”  (J=3.05).  Density  distribution 
(white:  60  kg  m-3  ;  black:  410  kg  m-3  ;  logarithmic  scale).  Sensor  B  is  placed  at  the 
end  of  the  inner  jet  potential  core  which  corresponds  to  different  location  in  cases  “N2” 
and  “N6”.  (b)  Case  “N8”  (J  =  9.3).  Density  distribution  (white:  60  kg  m-3  ;  black: 
maximum  ;  logarithmic  scale).  Sensor  C  is  located  in  the  recirculation  region. 
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5.3.1  Cases  “N2”  and  “N6” 


To  examine  the  spectral  information  it  is  natural  to  define  a  set  of  Strouhal  numbers  as¬ 
sociated  with  characteristic  scales  and  velocities  in  the  region  of  interest.  The  meaningful 
length  scales  in  this  study  are  the  inner  injector  diameter  di7  the  inner  lip  thickness  li  and 
the  annular  jet  thickness  he.  Two  velocities  characterize  the  system  :  the  inner  and  outer 
jet  injection  velocities  (ut  and  ue,  respectively).  The  resulting  Strouhal  numbers  and  the 
probes  used  to  measure  frequencies  are  shown  in  Tab.  4.  Probes  A  and  B  are  used  to 
evaluate  the  fundamental  frequencies  of  the  outer  and  inner  streams.  Probe  C  provides 
the  dominant  frequency  associated  with  eddy  structures  released  behind  the  injector  lip 
and  associated  with  the  outer  jet.  From  instantaneous  snapshots  (see  Fig.  12  for  exam¬ 
ple)  and  animations,  these  structures  are  found  to  originate  from  the  lip  corner  on  the 
outer-jet  side.  Positions  of  the  different  probes  are  shown  in  Fig.  12  (a).  Note  that  the 
axial  location  of  probe  B  is  not  the  same  in  cases  “N2”  and  “N6”,  because  this  probe 
is  placed  at  the  end  of  the  potential  core.  Power  spectral  densities  are  estimated  using 
Welch’s  method  of  periodograms.  Statistical  averaging  is  obtained  by  segmenting  the 
finite  data  set  in  M  overlapping  blocks  (a  50  %  overlap  is  used  to  augment  the  number 
of  periodograms  used  in  the  averaging  procedure).  Because  the  number  of  samples  is 
limited,  the  averaging  is  carried  out  over  M  =  8  blocks  in  the  outer  flow  and  on  M  =  4 
blocks  in  the  inner  jet  (power  is  concentrated  in  the  low  frequency  range  and  it  is  impor¬ 
tant  to  augment  the  resolution  to  distinguish  the  maximum  frequency).  The  frequency 
resolution  corresponding  to  the  different  sensors  is  given  in  Tab.  4. 


St  =  ( fl)/u 

l 

u 

Probe  (Fig.  12(a)) 

A/  [Hz] 

Ste 

he  (0.41  mm) 

Ue 

A 

N2:  375  ;  N6:  250 

Sf 

di  (0.51  mm) 

Ui 

B 

N2:  187  ;  N6:  125 

St1 

li  (0.54  mm) 

Ue 

C 

N2:  375  ;  N6:  250 

Table  4:  Definition  of  characteristic  Strouhal  numbers.  The  characteristic  length  and 
velocity  /  and  u  are  given  in  the  second  and  third  columns.  The  frequency  /  corresponds  to 
the  power  spectral  density  maximum  estimated  from  the  simulations.  A /  is  the  frequency 
resolution. 
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The  outer  and  central  streams  are  destabilized  at  a  few  diameters  ci,;  from  the  injector 
exit.  The  measured  frequency  for  the  outer  jet  (probe  A)  is  3  600  Hz  for  case  “N2”  and 
8  550  Hz  for  case  “N6”  (Fig.  13),  which  correspond  to  St^r2=0.25  and  St^r6=0.25.  The 
central  high-density  jet  presents  a  fundamental  instability  at  1  620  Hz  and  1  740  Hz  for 
cases  “N2”  and  “N6”  (Fig.  14)  respectively  (probe  B).  The  associated  Strouhal  numbers 
are  then  StjV2=0.34  and  Stjy6=0.26.  For  both  cases,  the  fundamental  mode  of  the  inner 
and  outer  streams  falls  in  the  typical  range  for  round  jets  (0.2  -  0.4). 


Outer  jet 


Figure  13:  Power  spectral  density  of  the  radial  velocity  in  the  outer  jet  (probe  A).  Aver¬ 
aging  is  carried  out  on  M  —  8  blocks  with  a  50%  overlap,  (a)  Cases  “N2”  (J=1.05),  (b) 
Case  “N6”  (J=3.05). 


Inner  jet 


Figure  14:  Power  spectral  density  of  radial  velocity  in  the  inner  jet  (probe  B).  Averaging 
is  carried  out  on  M  =  4  blocks  with  a  50%  overlap,  (a)  Cases  “N2”  (J=1.05),  (b)  Case 
“N6”  ( <7=3. 05). 


The  mixing  layer  between  the  two  jets  features  eddy  structures  generated  downstream 


32 


of  the  central  lip.  Their  dynamics  are  similar  to  those  observed  for  vortices  shed  be¬ 
hind  a  backward-facing  step.  The  characteristic  frequencies  are  found  to  be  1  620  Hz 
and  3600  Hz,  for  cases  “N2”  and  “N6”  (Fig.  15),  respectively.  Such  frequencies  corre¬ 
spond  to  Strouhal  numbers  St^— 0.15  and  St]y6=0.23  for  cases  “N2”  and  “N6,”  which 
are  slightly  higher  than  the  commonly  determined  values  for  backward-facing  step  in¬ 
stabilities  (St=O(0.1)  [70]).  This  could  be  attributed  to  a  reduction  of  the  “effective” 
lip  thickness  as  shear  is  increased.  It  turns  out  that  the  backward-facing  step  vortical 
structures  induce  a  back-flow  behind  the  lip  (Fig.  9).  High  density  fluid  from  the  central 
jet  is  entrained  and  recirculates  behind  the  lip.  One  may  then  consider  that  this  results  in 
a  reduced  “effective”  lip  thickness  which  intervenes  in  the  formation  of  the  vortical  struc¬ 
tures,  thus  limiting  their  spatial  development  and  increasing  the  instability  frequency. 
All  the  Strouhal  numbers  obtained  are  gathered  in  Tab.  5.  The  frequency  obtained  for 
case  “N2”  is  interestingly  close  to  that  associated  with  the  inner  jet  preferred  mode.  This 
observation  may  explain  why  large  scale  structures  are  exclusively  observed  in  this  case 

(Fig.  6). 


Inner  lip 


Figure  15:  Power  spectral  density  of  radial  velocity  in  the  mixing  layer  behind  the  inner 
lip  (probe  C).  Averaging  is  carried  out  over  M  =  8  blocks  with  a  50%  overlap,  (a)  Cases 
“N2”  ( J=1.05),  (b)  Case  “N6”  (J=3.05). 


5.3.2  Case  “N8”:  recirculation  regime 

The  mechanism  leading  to  mixing  of  the  high  density  jet  with  the  outer  stream  is  some¬ 
what  different  for  case  “N8”.  It  is  already  known  that  the  dense  jet  ends  abruptly  and 
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Case 

Ste  (probe  A) 

St1  (probe  B) 

St1  (probe  C) 

N2 

0.25 

0.34 

0.15 

N6 

0.25 

0.26 

0.23 

Tabic  5:  Characteristic  Strouhal  numbers  obtained  from  the  simulations  for  cases  “N2” 
and  “N6”. 

that  a  recirculation  bubble  is  established  beyond  that  point.  As  done  previously  for  cases 
“N2”  and  “N6”,  characteristic  length  scales  used  to  define  the  relevant  Strouhal  numbers 
are  given  in  Tab.  6,  and  the  probes  used  to  characterize  this  case  are  shown  in  Fig. 
12(b).  Probe  C  is  inside  the  recirculation  region  and  is  located  at  x  =  xu.  This  case  is 
highly  turbulent  and  features  a  broadband  spectrum  so  that  temporal  data  are  arranged 
in  M  =  16  blocks  to  ensure  a  reasonable  statistical  convergence.  A  50  %  overlap  is  used 
to  obtain  the  maximum  possible  averaging  number.  The  resulting  frequency  resolution 
is  in  this  case  1  000  Hz.  This  value  may  seem  quite  high,  but  should  be  compared  to  the 
spectral  range  which  extends  to  25  000  Hz. 


St  =  ( fl)/u 
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u 

probe  (Fig.  12(b)) 

00 
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he  (0.41  mm) 

ue 
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JiN  8 

(0.54  mm) 

ue 
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Q+b 
oiN  8 

di  +  h  (0.51+0.54  mm) 

Ue 

C 

Table  6:  Strouhal  numbers  definitions.  I  and  u  are  the  characteristic  length  and  velocity 
used  to  define  the  Strouhal  numbers.  /  represents  the  frequency  obtained  from  the 
simulations  and  A /  is  the  frequency  resolution. 

The  dominant  frequency  for  the  outer  jet  (probe  A)  is  21  600  Hz  (Fig.  16),  which 
corresponds  to  St^8=0.25,  as  also  obtained  for  the  two  other  cases.  The  inner  jet  is 
notably  affected  by  the  recirculation  bubble.  Its  longitudinal  development  is  limited  and 
mass  is  ejected  in  the  radial  direction  in  the  form  of  thin  layers  appearing  in  Fig.  12(b). 
The  resulting  structure  is  similar  to  a  counterflow.  This  motion  leads  to  the  formation 
of  a  recirculation  ring  of  high  density  fluid  behind  the  inner  lip  with  a  characteristic 
size  of  the  order  of  the  lip  thickness.  Vortex  roll-up  is  observed  behind  the  inner  lip, 
in  the  mixing  layer  established  between  the  inner  and  outer  jets.  These  structures  are 
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similar  to  those  observed  in  cases  “N2”  and  “N6”,  but  they  strongly  interact  with  higher 
density  layers  emerging  from  the  central  jet,  and  their  spatial  development  is  limited  by 
the  recirculation  region.  The  corresponding  Strouhal  number  is  StlN8= 0.26  (Fig.  17(a)). 
These  eddy  structures  lead  to  the  formation  of  the  recirculation  region.  The  measured 
frequency  corresponds  to  StbN8= 0.18  (Fig.  17(b)).  The  dynamics  of  this  bubble  are 
similar  to  those  observed  behind  a  bluff  body  terminated  by  a  flat  face,  for  which  the 
characteristic  Strouhal  number  is  close  to  0.2.  The  Strouhal  numbers  obtained  in  this 
analysis  are  gathered  in  Tab.  7. 

Outer  jet 


Frequency  [Hz] 

Figure  16:  Case  “N8”  (J=9.3).  Power  spectral  density  of  pressure  on  Probe  A.  Averaging 
is  carried  out  over  M  =  16  blocks  with  50%  overlap. 


Probe  B  Recirculation  bubble 


Figure  17:  Case  “N8”  (J=9.3).  (a)  Power  spectral  density  of  pressure  on  Probe  B,  (b) 
Power  spectral  density  of  pressure  in  the  recirculation  bubble.  Averaging  is  carried  out 
over  M  =  16  blocks  with  50%  overlap. 
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SteN g  (probe  A)  StlN 8  (probe  B) 

StbN 8  (probe  C) 

0.25  0.26 

0.18 

Table  7:  Characteristic  Strouhal  numbers  obtained  from  the  simulations  for  case  “N8”. 

6  Effects  of  a  transverse  acoustic  modulation 

Effects  of  an  acoustic  modulation  are  now  investigated  in  case  “N6”.  This  operating 
point  is  chosen  because  the  core  is  sufficiently  long  to  allow  an  examination  of  external 
modulation  effects.  The  methodology  used  to  impose  an  acoustic  excitation  is  described 
in  Sec.  6.1.  Results  are  then  compared  with  experimental  observations  in  Sec.  6.2  and 
discussed  in  Sec.  6.3. 

6.1  Modulation  methodology 

An  acoustic  modulation  is  obtained  by  imposing  harmonic  modulations  of  the  normal 
velocity  at  the  outer  boundaries  using  a  method  devised  by  Rey  et  al.  [52]  (Fig.  18). 
The  modulation  has  the  form  : 

=  sin(27r/t)  v2  =  sin(27r/t  +  0)  (18) 

where  V\  and  v2  are  the  normal  velocities  at  the  lateral  boundaries,  /  is  the  modulation 
frequency  set  equal  to  3  000  Hz,  0  designates  the  phase  between  the  imposed  modulations. 
All  other  boundaries  are  treated  as  in  the  previous  cases  (in  the  absence  of  modulation). 
This  technique  is  first  validated  by  simulating  an  acoustic  perturbation  in  the  absence  of 
the  central  flow.  Longitudinal  distributions  of  the  rms  pressure  fluctuations  are  shown  in 
Fig.  19.  When  0  =  0  (Fig.  19(a)),  the  normal  velocities  at  the  boundaries  are  in  phase, 
and  this  yields  a  maximum  pressure  fluctuation  on  the  injector  axis.  A  pressure  probe 
placed  at  5  mm  from  the  injector  exit  confirms  this  behavior.  As  expected,  pressure 
fluctuations  are  observed  with  a  2  %  peak-to-peak  oscillation  with  respect  to  the  mean 
chamber  pressure  while  the  transverse  velocity  remains  unperturbed.  When  the  normal 
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velocities  at  the  boundaries  are  in  phase  opposition  (0  =  7 r),  the  pressure  vanishes  on 
the  centerline  and  the  transverse  velocity  oscillation  reaches  its  maximum  (Fig.  19(b)). 
Spectral  analysis  also  indicates  that  the  perturbed  held  is  dominated  by  the  imposed 
modulation  frequency  and  that  the  harmonic  content  is  negligible. 


Figure  18:  Longitudinal  slice  in  the  computational  domain  showing  the  methodology 
retained  to  impose  acoustic  perturbation. 


Figure  19:  Results  obtained  under  acoustic  modulation  in  the  absence  of  how.  (a)  RMS 
pressure  fluctuations  for  acoustic  perturbations  in  phase  (0  =  0).  (b)  RMS  pressure 
huctuations  for  acoustic  perturbations  in  phase  opposition  (0  =  7r). 


The  simulated  cases  are  defined  in  Tab.  8.  Injection  conditions  correspond  to  case  “N6”. 
Modulation  amplitudes  of  1  and  2%  of  the  mean  chamber  pressure  are  imposed.  These 
two  cases  are  studied  for  0  =  0  and  0  =  7 r,  respectively,  corresponding  to  a  pressure 
antinode  and  a  pressure  node  on  the  central  axis. 
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Case 

N6p0 

N6^ 

N6}*, 

N6^ 

Ppp/Po 

0 

2.0  % 

0 

2.0  % 

7 T 

1.0  % 

0 

1.0  % 

7T 

Table  8:  Test  cases  with  acoustic  modulation. 

6.2  Calculated  flow  dynamics  and  experimental  data 

When  acoustic  modulations  are  in  phase,  pressure  fluctuations  are  produced  at  the  in¬ 
jector  outlet  at  the  modulation  frequency.  This  leads  to  oscillations  of  the  axial  velocity 
at  the  imposed  frequency.  Such  variations  in  the  longitudinal  velocity  induce  a  periodic 
shedding  of  ring- like  patterns  (Fig.  20).  Structures  are  created  behind  the  inner  lip  and  in 
the  outer  mixing  layer.  These  structures  develop  in  an  axi-symmetric  fashion  producing 
annular  vortices.  The  vortices  produced  in  the  inner  jet  lip  entrain  fluid  from  the  inner 
jet,  involving  mass  from  the  inner  jet  in  the  vortex  roll-up.  These  results  are  qualitatively 
consistent  with  experimental  observations  shown  in  Fig.  22(a),  where  a  varicose  motion 
is  clearly  identifiable.  When  acoustic  perturbations  are  in  phase-opposition  (Fig.  21), 

K  ¥  -  if-  - 

□  □  □ 


Figure  20:  Case  “N6^0”  (J=3.05).  Longitudinal  slices  of  density  (white:  60  kg  m-3  ; 
black:  maximum  ;  logarithmic  scale).  These  images  correspond  to  three  successive  times 
in  one  period. 


the  resulting  transverse  velocity  perturbation  induces  a  sinusoidal  motion,  as  observed 
experimentally  in  Fig.  22(b);  and  vortices  are  shed  in  an  asymmetric  fashion. 

The  centerline  density  and  axial  velocity  are  plotted  in  Fig.  23.  Acoustic  perturba¬ 
tions  significantly  reduce  the  dense  core  length  and  increase  the  amplitude  of  the  low 
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Figure  21:  Case  “N (J=3.05).  Longitudinal  slices  of  density  (white:  60  kg  m~3  ; 
black:  maximum  ;  logarithmic  scale).  These  images  correspond  to  three  successive  times 
in  one  period. 


Figure  22:  Experimental  shadowgraphs  for  case  “N6”.  (a)  Acoustic  modulation  is  in 
phase,  (b)  Acoustic  modulation  is  in  phase-opposition. 
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velocity  region  just  behind  the  central  jet.  Characteristic  length  ratios  between  acous¬ 
tically  perturbed  cases  and  the  reference  case  “N6”  are  given  in  Tab.  9.  These  ratios 
are  in  agreement  with  experimental  measurements.  It  is  interesting  to  note  that  the 
cases  corresponding  to  phase  opposition  produce  the  shortest  dense  core.  It  is  also  found 
that  multiplication  of  the  perturbation  amplitude  by  a  factor  of  two  does  not  lead  to  a 
proportional  reduction  of  the  dense  core. 


Figure  23:  (a)  Density  profile  on  the  centerline,  (b)  Axial  velocity  profile  on  the  centerline. 


Case 

N6^0 

N6f° 

N6}i 

N6i\ 

0— 7T 

70.5  / 70.5 

L  p  /  L  p,noac 

0.60 

0.34 

0.70 

0.45 

p.i/p.i 
v p  /  vp,noac 

0.53 

0.36 

0.66 

0.45 

Experiments 

- 

- 

0.68 

0.55 

Table  9:  Influence  of  acoustic  modulations  on  the  characteristic  length  in  case  “N6”. 


6.3  Spectral  analysis  of  acoustic  modulation  effects 

The  response  of  the  flow  to  the  external  modulation  may  now  be  discussed  by  examining 
the  spectral  content  of  velocity  fluctuations.  This  is  done  by  plotting  spectral  maps  of 
these  fluctuations  along  three  axial  lines  in  the  flow.  The  first  (line  C)  corresponds  to 
the  centerline,  the  second  (line  B)  is  centered  on  the  lip,  the  third  (line  A)  is  centered  on 
the  outer  jet  annulus  (see  Fig.  24). 
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Axis  A 


Axis  B 
Axis  C 


Figure  24:  Case  “N6|^0”  ( <7=3. 05).  Positions  of  the  three  axial  lines  used  in  the  spatial 
spectral  analysis. 

6.3.1  Case  “N6JJ,” 

In-phase  excitation  induces  axial  velocity  modulations  on  the  centerline.  This  is  clearly 
seen  in  Fig.  25(a),  26(a)  and  27(a),  where  the  3  000  Hz  fluctuations  are  well  retrieved. 
The  amplitude  of  the  axial  velocity  fluctuation  is  nevertheless  very  low  in  the  high-density 
region  (Fig.  28),  in  the  potential  core  of  the  inner  jet  (;x  <1.05  dj,  defined  by  the  region 
along  the  centerline  axis  where  p  >  0.99 p,,  see  Fig.  23(a)).  This  region  of  the  flow  is  less 
affected  by  the  axial  modulation  than  the  surround  streams.  The  response  of  the  jets  to 
this  modulation  is  determined  by  examining  the  transverse  velocity.  Both  the  outer  jet 
and  the  flow  behind  the  inner  lip  show  a  strong  response  at  3  000  Hz  (Fig.  25(b)  and 
26(b)).  It  means  that  the  axial  perturbation  promotes  a  destabilization  of  these  jets  and 
vortex  roll-up  at  the  same  frequency.  This  is  not  the  case  along  the  inner  jet  centerline, 
where  the  transverse  response  at  3  000  Hz  is  not  distinguishable  (Fig.  27(b)).  Since  the 
geometry  and  the  perturbation  are  both  axisymmetric  around  this  axis,  such  a  velocity 
motion  is  naturally  limited.  The  dominant  frequency  after  the  potential  core  ( x  >1.05  d,) 
is  close  to  1  000  Hz.  The  spectral  content  of  the  inner  jet  is  shown  in  Fig.  28  at  x  =  dj, 
at  the  end  of  the  potential  core.  The  axial  velocity  is  dominated  by  the  modulation  at  3 
000  Hz,  with  a  low  amplitude,  but  also  features  a  lower  frequency  component  around  1 
200  Hz.  In  contrast  with  the  outer  jet,  the  radial  velocity  responds  at  half  of  the  imposed 
frequency  (around  1  500  Hz).  These  fluctuations  are  nevertheless  small  compared  to  those 
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measured  in  the  gaseous  part  of  the  flow.  As  a  consequence,  the  high  density  part  of  the 
inner  jet  seems  to  be  virtually  insensitive  to  the  acoustic  modulation.  The  frequency  that 


dominates  at  the  end  of  the  potential  core  is  close  to  half  of  the  modulation  frequency. 
This  frequency  is  also  observed  in  the  analysis  of  density  fluctuations. 
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Figure  25:  Case  “N6^0”  (J=3.05).  Spectral  map  along  line  A  (r  =  di/2  + +  he/2).  (a) 
Axial  velocity,  (b)  Radial  velocity. 
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Figure  26:  Case  “N6|^0”  (J=  3.05).  Spectral  map  along  line  B  (r  =  di/2  +  k/2).  (a) 
Axial  velocity  and  (b)  Radial  velocity. 


6.3.2  Case 

A  transverse  velocity  modulation  is  again  applied  at  3  000  Hz,  but  the  phase  is  now  equal 
to  7 r.  As  in  the  previous  case  “N6|^0”,  this  transverse  modulation  affects  all  the  jets  as 
shown  in  Fig.  29  (b),  30  (b)  and  31  (b).  The  amplitude  of  this  transverse  perturbation 
is  much  lower  in  the  high  density  part  of  the  inner  jet  (Fig.  32  and  31(b)  for  x  <  d/). 
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Figure  27:  Case  “N6|=0”  (J—  3.05).  Spectral  map  along  line  C  (r  =  0).  (a)  Axial 
velocity,  (b)  Radial  velocity.  The  dashed  line  represents  the  axial  position  at  which  a 
power  spectral  density  is  shown  in  Fig.  28. 


Figure  28:  Case  “N6^0”  ( J —  3.05).  Power  spectral  density  of  axial  ( — )  and  radial  (— 
— )  velocity  in  the  inner  jet  centerline  at  x  =  di  (end  of  the  potential  core). 
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Axial  velocity  perturbations  are  induced  in  the  outer  jet  and  behind  the  inner  lip  (Fig. 
29(a)  and  30(a)),  associated  with  vortex  roll-up  at  the  same  frequency.  Again,  the  inner 
jet  is  essentially  unaffected,  since  the  measured  amplitudes  are  much  lower  than  for  the 
two  other  cases  (Fig.  31(a)).  This  stream  is  dominated  by  an  intermediate  frequency.  A 
component  at  1  500  Hz  is  observed  just  at  the  inner  injector  exit.  This  is  confirmed  in  Fig. 
32,  where  the  transverse  modulation  at  3  000  Hz  is  visible,  but  leads  to  a  longitudinal 
response  of  the  jet  close  to  1  500  Hz.  Even  though  this  frequency  is  clearly  identifiable 
in  Fig.  31(a)  and  32,  its  amplitude  is  lower  than  that  measured  behind  the  lip  or  in  the 
outer  jet.  As  previously  observed  for  in-phase  modulation,  the  sensitivity  of  the  inner 
jet  to  the  external  modulation  is  limited.  Again,  the  frequency  which  dominates  at  the 
end  of  the  potential  core  {x  <0.55  dj,  see  Fig.  23(a))  is  close  to  half  of  the  modulation 
frequency  and  is  also  observed  in  density  fluctuations. 

Axis  A 


Axial  velocity  Radial  velocity 


Figure  29:  Case  “NO^.”  (J=  3.05).  Spectral  map  along  line  A  (r  =  di/2  +  k  +  he/2). 
(a)  Axial  velocity  and  (b)  Radial  velocity. 


6.4  Discussion 

It  is  found  that  the  inner  jet,  when  submitted  to  an  acoustic  modulation  at  3  000  Hz, 
exhibits  fluctuations  at  the  subharmonic  frequency  of  1  500  Hz.  This  may  be  explained 
by  considering  the  stability  properties  of  the  central  jet  and  by  analyzing  the  acoustic 
transmission  from  the  outer  medium  into  the  high  density  region. 
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Figure  30:  Case  “NO^y’  (J=  3.05).  Spectral  map  along  line  B  (r  =  di/2  +  U/2).  (a) 
Axial  velocity  and  (b)  Radial  velocity. 
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Figure  31:  Case  “N6^2^.”  (J—  3.05).  Spectral  map  along  line  C  (r  =  0).  (a)  Axial 
velocity  and  (b)  Radial  velocity.  The  dashed  line  represents  the  axial  position  at  which 
a  power  spectral  density  is  shown  in  Fig.  32. 


Figure  32:  Case  “NO2)^’  (J=  3.05).  Power  spectral  density  of  axial  ( — )  and  radial  (— 
—)  velocity  on  the  inner  jet  centerline  at  x  —  0.5 A  (end  of  the  potential  core). 
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It  is  first  interesting  to  situate  the  modulation  frequency  with  respect  to  the  range  of 
unstable  frequencies  characterizing  the  jet  and,  more  specifically,  with  respect  to  the 
preferred  frequency  of  the  inner  jet.  It  is  known  (see  for  example  Birbaud  et  ah  [3]  and 
references  in  this  article)  that  when  the  modulation  frequency  is  well  above  the  preferred 
mode  frequency,  the  flow  is  not  affected  by  an  external  acoustic  modulation.  This  is 
best  expressed  in  terms  of  a  Strouhal  number  which  can  be  calculated  by  considering 
the  acoustic  modulation  frequency,  the  characteristic  length  and  velocity  associated  with 
each  jet.  The  resulting  “acoustic”  Strouhal  numbers  are  given  in  Tab.  10.  While  St®c  and 
Stlac  are  lower  or  close  to  the  natural  Strouhal  number  of  the  flow  for  case  “N6”,  Sttc  is 
higher  than  the  one  measured  previously  without  modulation.  As  a  consequence,  the  jet 
does  not  respond  to  the  modulation  frequency,  but  oscillates  at  a  subharmonic  frequency 
(1  500Hz)  which  is  close  to  the  preferred  frequency  measured  previously  in  the  absence 
of  acoustic  modulation  (Fig.  14). 


St  =  ( U)/u 

length  scale  l 

charac.  vcl.  u 

value 

steac 

he  (0.41  mm) 

ue 

0.087 

stlac 

(0.54  mm) 

ue 

0.115 

Stlac 

di  (0.51  mm) 

Ui 

0.45 

Table  10:  Definition  of  the  Strouhal  numbers  used  in  this  study.  /  represents  the  modu¬ 
lation  frequency  (3  000  Hz). 

Another  possible  explanation  of  the  absence  of  modulation  at  the  imposed  frequency  may 
be  obtained  by  considering  the  change  in  acoustic  impedance  between  the  outer  stream 
and  the  inner  jet.  It  is  known  that  a  large  impedance  mismatch  reduces  the  transmission 
of  waves  across  the  interface.  Acoustic  impedances  gathered  in  Tab.  11  indicate  that  a 
large  variation  takes  place  at  the  inner  jet  boundary.  Under  these  circumstances  acous¬ 
tic  wave  transmission  into  the  jet  is  reduced,  explaining  the  absence  of  the  modulation 
imposed  to  the  system  from  the  outer  region. 

The  previous  arguments  provide  some  clues  on  the  observed  response  of  the  coaxial 
flow  to  an  externally  imposed  transverse  acoustic  motion.  The  strong  reduction  of  the 
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Conditions 

Impedance  ( Z  =  pc) 

Chamber 

17  701 

Outer  jet 

19  663 

Inner  jet 

98  700 

Table  11:  Specific  acoustic  impedance  in  the  chamber,  inner  and  outer  jets. 

inner  jet  length  is  certainly  linked  to  the  strong  response  of  the  outer  jet  to  the  acoustic 
modulation,  which  is  not  observed  for  the  inner  stream.  The  large  vortices  generated 
in  the  annular  jet  in  turn  enhance  mixing  between  the  outer  and  inner  jets,  and  thus 
enhance  mixing  of  the  inner  fluid  with  its  surroundings,  leading  to  a  reduction  of  the 
dense  core  length. 

7  Conclusion 

The  structure  of  non-reactive  transcritical  and  supercritical  jets  is  investigated  by  exam¬ 
ining  the  flow  formed  by  a  coaxial  injector  operating  at  high  pressure.  This  unit  injects 
low  temperature  dense  nitrogen  surrounded  by  an  outer  stream  of  higher  temperature 
nitrogen.  The  analysis  relies  on  large  eddy  simulations  which  are  compared  with  ex¬ 
perimental  results.  The  natural  development  of  the  coaxial  flow  is  examined  for  a  set 
of  momentum  flux  ratios.  It  is  shown  that  the  flow  patterns  observed  experimentally 
are  suitably  retrieved  and  that  a  good  correspondence  between  experimental  results  and 
simulations  is  obtained.  It  is  shown  in  particular  that  when  the  momentum  flux  ratio 
exceeds  a  critical  value,  the  central  jet  is  terminated  in  an  abrupt  fashion  at  a  short 
distance  from  the  injector  exhaust.  Calculations  indicate  that  this  is  linked  to  the  forma¬ 
tion  of  a  recirculation  region  located  at  a  short  distance  from  the  injection  section.  This 
study  is  complemented  with  a  spectral  analysis  of  the  velocity  components  in  the  injector 
nearfield.  It  is  shown  that  the  frequency  content  is  governed  by  Strouhal  numbers  which 
follow  standard  rules  for  jets  and  wakes.  The  effect  of  an  external  acoustic  modulation  is 
then  investigated  by  imposing  in-phase  and  out-of-phase  perturbations  at  the  boundaries 
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of  the  domain.  This  is  used  to  study  the  effect  of  transverse  waves  on  the  dark  core 
in  the  jet.  A  strong  reduction  of  the  jet  length  is  clearly  observed,  both  in  experiment 
and  numerics.  The  spectral  content  of  fluctuations  in  the  coaxial  streams  is  analyzed. 
In  the  case  considered  which  corresponds  to  a  Strouhal  number  (based  on  the  inner  jet 
velocity  and  diameter)  about  equal  to  twice  the  Strouhal  number  corresponding  to  the 
preferred  mode  of  the  jet,  it  is  found  that  the  inner  jet  responds  with  a  limited  amplitude 
to  the  external  modulation  but  features  a  subharmonic  motion  at  half  of  the  modulation 
frequency.  In  contrast,  the  gaseous  part  of  the  flow  is  strongly  perturbed  by  the  acoustic 
modulation.  Large  scale  vortices  are  released  from  the  outer  annulus  at  the  modulation 
frequency.  It  is  suggested  that  the  relative  insensitivity  of  the  inner  flow  is  linked  to  the 
fact  that  the  frequency  is  above  the  range  of  unstable  frequencies  of  the  central  jet.  It  is 
also  found  that  the  specific  acoustic  impedance  of  the  central  jet  is  much  higher  than  that 
corresponding  to  the  surrounding  stream  and  that  this  may  explain  why  the  modulation 
is  not  transmitted  to  the  central  core. 
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A  Injection  conditions  for  the  computed  cases 


Outer  jet  injection  conditions 


Case 

Te  [K] 

me  [g.s-1 

pe  [kg.m  3[ 

ue 

m.s  1 

Ree 

N2 

152 

1.57 

101 

5.95 

41  600 

N6 

183 

2.69 

73.1 

14.1 

64  800 

N8 

192 

6.46 

68.3 

36.2 

151  000 

Table  12:  Injection  conditions  corresponding  to  the  cases  “N2”,  “N6”  and  “N8”  for  the 
outer  jet.  m  denotes  the  mass  flow  rate. 


Inner  jet  injection  conditions 


Case 

Ti  [K] 

m 

g.s  1 

Pi  [kg.m  3[ 

Hi 

m.s  1 

Rei 

N2 

117 

0.289 

590 

2.4 

15  000 

N6 

126 

0.292 

420 

3.4 

25  000 

N8 

128 

0.295 

220 

6.6 

52  000 

Table  13:  Injection  conditions  corresponding  to  the  cases  “N2”,  “N6”  and  “N8”  for  the 
inner  jet.  m  denotes  the  mass  flow  rate. 
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